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One of the broadest concepts of measurement in quantum theory is the generalized measurement. 
Another paradigm of measurement — arising naturally in quantum optics, among other fields — is that 
of continuous-time measurements, which can be seen as the limit of a consecutive sequence of weak 
measurements. They are naturally described in terms of stochastic processes, or time-dependent 
random variables. We show that any generalized measurement can be decomposed as a sequence 
of weak measurements with a mathematical limit as a continuous stochastic process. We give an 
explicit construction for any generalized measurement, and prove that the resulting continuous 
evolution, in the long-time limit, collapses the state of the quantum system to one of the final states 
generated by the generalized measurement, being decomposed, with the correct probabilities. A 
prominent feature of the construction is the presence of a feedback mechanism — the instantaneous 
choice weak measurement at a given time depends on the outcomes of earlier measurements. For 
a generalized measurement with n outcomes, this information is captured by a real n-vector on an 
n-simplex, which obeys a simple classical stochastic evolution. 



I. INTRODUCTION 

The first and simplest definition of measurement that one 
learns in a first course on quantum mechanics involves pro- 
jection operators which sum up to the identity [lj. Usually, 
these projective operators come from a resolution of the iden- 
tity generated by a particular observable, represented by a 
Hermitian operator. A generalization is needed when com- 
posite systems are considered — a projective measurement on 
the whole system is not usually a projective measurement on 
a part of the system. Every generalized measurement can be 
performed by adding an ancilla to our given system, doing 
an appropriate unitary transforming both system and ancilla, 
followed by a projective measurement on the ancilla The 
dimension of the ancilla must be at least equal to the num- 
ber of outcomes from the measurement and the number of 
orthogonal projectors on the ancilla. 

A projective measurement can be carried out as a sequence 
of weak generalized measurements Here, weak means that 
after each step in the sequence, the system is disturbed by a 
small amount, but yields only a small amount of information. 
Each step also takes a certain finite amount of time, and thus 
so will the final strong projective measurement. This pic- 
ture of a measurement, as a continuous sequence of infinites- 
imal steps taking a finite amount of time, contrasts with the 
usual assumptions we make about abstract measurements — 
that they are instantaneous and thus strong. 

While this construction works for projective measurements, 
it cannot be immediately adapted to strong generalized mea- 
surements. We wish to construct a measurement procedure 
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that is continuous in time, but that at long times produces 
the same result as a given strong measurement: that is, in 
the end the system is in one of the correct final states with 
the correct probabilities. We will prove that such a contin- 
uous decomposition exists for any generalized measurement, 
and give an explicit construction. Such a continuous proce- 
dure can serve a variety of purposes. For certain problems, a 
continuous description is useful: we are able to take deriva- 
tives and use other analytical tools of calculus [H, 0], apply 
methods from quantum filtering theory and quantum feed- 
back control. On the level of physical reality, it is difficult to 
argue whether one particular description of a measurement is 
more fundamental (if any of them is). For much of the time 
since the discovery of quantum mechanics it made sense to 
treat measurements as instantaneous because those were the 
kind of measurements we were able to perform. With more 
recent advances in quantum optics and atomic physics, where 
quantum systems can be continuously monitored, new ways of 
describing the evolution of the system had to be developed — 
models such as quantum trajectories or decoherent histories 
H 1, 0, H, II E3, El El- That there are numerous ways to 
think about measurements in quantum mechanics, which are 
not necessarily contradictory, but on the contrary are in some 
sense complementary. 

The purpose of this paper is to construct a quantum 
continuous stochastic process — a family of consecutive weak 
measurements — which governs the evolution of the state of 
the quantum system, and for which the final result (at long 
times) is the same as a that of a specified (strong) general- 
ized measurement. We will give an explicit construction of 
such a process, and prove that it does indeed have the correct 
long-time behavior. As we will see, in general this requires 
that the choice of measurements at later times will depend 
on the measurement results at earlier time, so this continuous 
measurement procedure must include feedback. 

In Section II we first consider the discrete case of succes- 
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sive weak but finite measurements as a good and fairly simple 
introduction to what happens in the continuous case. After 
that, in Section III we derive a stochastic differential equation 
for the continuous stochastic process implementing a projec- 
tive measurement. This construction is simplified by the fact 
that the measurement operators all commute, and the final 
states are all orthogonal, but it will show the approach to be 
followed in the more general case. We prove that our stochas- 
tic process really does generate the results of the generalized 
measurement as time goes to infinity. Finally, in Section IV 
we show how the same thing can be done for any generalized 
measurement. This does not require that the measurement 
operators commute, nor that the final states be orthogonal, 
and in general requires feedback from the results of earlier 
measurements to the choice of later measurements. We give 
an example, in the case when the measurement operators do 
commute, of how generated quantum state diffusion equations 
can be derived from our equations. Section V briefly summa- 
rizes our results and conclusions. 



II. REPEATED WEAK MEASUREMENTS FOR 
PROJECTIVE MEASUREMENTS 

Let's consider a projective measurement on the system. 
The measurement operators are denoted by Pi and they have 
the following properties 



p. p. _ X. . p. 



(1) 
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We would also assume that initially our system is in state 
\ipo) ■ After we perform the projective measurement ([1]) on 
the system, it states collapses to one of the following states 



Pi \4>o) 



(2) 



with probability p® — (ipo\ Pi \ipo) respectively. 

We now give the discrete procedure which decomposes the 
projective measurement into a sequence of weak measure- 
ments. We denote by A n the classical state space 



A n = jxe (o,i)"cl" 

and by A" - the closure of A" 

A n = \ x S [0, 1]™ C R n 



Let {x(k) ^ A n > k — 1, .., n} be n fixed points in A" with the 
property 



E 



X( k ) 



(3) 



where e = (i, -) 6 A n . These n points x^ will repre- 
sent the outcomes of the weak measurements that we do, and 



we will refer to them as the fundamental steps for our dis- 
crete stochastic process. The measurement operators that we 
construct out of them are 



(fc) 



(4) 



with the completeness condition satisfied: 

n n 

EVw= E • 

k — l fc,i,j — 1 



= E \J x \k) x U^j^ = E x hA 

k,i — l k,i — l 

n 

= n£VA = J. 

i-1 

This means that the measurement operators §4^ define a com- 
plete quantum measurement. If we measure the system using 
the above measurement operators over and over again, after a 
long enough time (strictly speaking in the limit of time going 
to infinity) the system is bound to end up in one of states ((2} 
with the right probability. After every measurement we get 
a measurement outcome (k), k G {1, n} so after s steps we 
have a sequence (fci, &2, k s ) of outcomes. The state of the 
system at that time has changed to 

= \i>(k s fc,.)) = . 1 JV (fcs) ..JV (fcl) |Vo) 



■v/Pfcs.-fci 



E 



2 n 



(5) 



where 



Pk 3 ... kl = {MN{ kl y..Nl ks) N (ks y..N (kl) \4, ) 

71 

= E^CM-^C*!) (V'ol-Pi IV'o) 

71 

Ei i 

X {ks)— X {hi)Pi ■ 

i=l 

Equation ([5]) can be rewritten as 

n 



with 



Because 



E^ = 



Pk s ...k 1 



ETi i i 

i=l X {k s )-- x {k 1 )Pi 



(6) 



(7) 



= 1 



it follows that x 3 £ A". In these new coordinates the initial 
state of the system is just |"0o) = X^"=i v^o 1$*) wnere x\, = 
p°. Another way of writing equation ((5]l is 



(8) 
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with 



En 
1=1 1 



(9) 



"(fcl) 



Here again x s £ A™. Then |t/jo} is given by ([8]) with x l — —. 

Let us stress here that both formulas © and ([8| imply that 
the state of the system changes after each weak measurement, 
generated by Q, and that there is a correspondence between 
the evolving vector x s = (x s , x") or x s — (x s , x% ), the 
time evolution of each given by formula and © respec- 
tively, and the evolving state of the system \ip s ) . 

Another observation is that instead of keeping track of the 
sequence of outcomes (fa, fa, ■•■) we can equally well use the 
sequence of points [xx,Xi, ...) or (xi,Xa,—) m A™. There is 
a one-to-one correspondence between sequences of outcomes 
and sequences of permissible points in A™. By a "permissible 
point" we mean a point for which there exists at least one 
sequence of outcomes (fa, fa, fa) that generates that point, 
by formula $7} or ([9]) respectively, given some initially chosen 
and fixed fundamental steps {£(&)}• To prove that all three 
type of sequences carry the same information we use the fact 
that A™ can be given a group structure: 

Definition 1 A™ is an abelian group with multiplication * : 
A" x A" — > A" where x-ky for x, y g A™ is defined as 



(x*yY 



(10) 



We also define the Hadamard product x o y of two vectors x 
and y of the same dimension by 

(xoy) = x y , 

and the trace Tr(x) of a vector x, which is defined as 

n 

Tr(x) = y]x i . 



Using these notations, equation (|10[) reads as 



x*y = 



xoy 
Tr(x o y) 



The identity e of the group A n is e = (1/n, 1/n), and the 
inverse a; -1 of an element x is 



Tr(x) 



where x ■ 



Now it is easy to see that given a sequence of permissible 
points (xi, X2, x s ), generated by the measurement proce- 
dure with some pre-fixed fundamental steps, every element of 
the sequence (xi,X2X^ 1 , x s xj_i) is one of the fundamental 
steps K(m for some k. Here x~ x is the inverse of x with respect 
to the group multiplication. So (x\, xix^ , x s k~2i) — 
(x( kl ),X( k2 ), ...,X( ks )) for some (fa,..., fa) chosen appropri- 
ately, which is exactly the sequence of outcomes. It's as easy 
to prove the same for a sequence (Si, X2, —,Xa)- 

The value of keeping track of sequences of permissible 
points is that we need only the last element x s of the sequence 
and the last outcome i 3 +i from the (s + l)-th measurement 
in order to find the state of the system \ip s +i) at time s + 1: 
it is given by formula |[S} with 



Xg-\-l — Xs * X^i 



Xs O X(j B+1 ) 

Tr(xs o a?(j , j)) 



(11) 



We can make use of this result by making a correspondence 
between the measurement procedure and a classical discrete 
stochastic process. Let (fi, J, P) be a probability space. We 
define the stochastic process as x : No X fl — > PA™ by giving 
its distribution law. Here PA™ is the set of permissible points 
for a concrete choice of fundamental steps Xn.) and No = 
{0, 1, 2, ...} is the time set. Usually we will denote x(s, u) as 
x s (u>) for u g fl. The probability distribution of x 3 thought 
of as random variable, is d s (x) — ^2p kl ,...,k s where the sum 
is over all sequences of outcomes (fa, fa) for which 



En 
1= 



(12) 



or in other words, for which the last entry in the sequence 
with s elements is exactly x. If there is no such sequence with 
s elements then d s (x) = 0. The reinterpretation of formula 
in these new terms is that the process x is Markovian. 



III. CONTINUOUS PROCESS FOR 
PROJECTIVE MEASUREMENTS 

We will now construct a continuous quantum measure- 
ments process which decomposes the strong projective mea- 
surement ([T}. Unlike the discrete case, we are now not re- 
stricted to performing the same measurement at every time 
step (in a sense, the fundamental steps can be chosen to be 
different). These measurements can be arbitrary, but they 
must depend smoothly on the state of the system at the cur- 
rent moment and not explicitly on the time. (This will be 
vital when we take the next step to decomposing generalized 
measurements, where the measurements must be different at 
different steps.) 

Let's make these ideas precise. As in the discrete case, we 
will concentrate on the evolution of the classical state ana- 
log of the system X t G A™ rather than its quantum state 
\ipt) ■ The measurement on the system is now given by a 
Lesbegue-measurable set D(x;8) C A™, of vectors which is a 
neighborhood of e £ A™, and satisfies a condition analogous 
to ©: 



/ 



zdht(z) 



dh*(z) = eA(x; 5) 



(13) 



D(x;S) 



with A(x; 8) being the (n — l)-dimensional volume of D(x; 5) 
and h+(x) the Haar measure on A™. The other condition 
needed is 



lim d(x: 5) 

6^0 



lim sup d(z, 

<5 -*° zGD(x;S) 







(14) 



where d(z, x) is the standard distance function on A™ thought 
of as a subset of R™. The set D(x; 5) depends on the current 
state of the system x and on a small positive parameter (0 < 
5 <g 1) which at some point we will let go to zero, thus getting 
a continuous stochastic process. The measurement operators 
are given as in Q: 

N(z, x; 5) = V2 Vz^Pi with z G D(x; 6) (15) 
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where the completeness condition is given by 



I 



D(x;S) 



n n 

/ 

<=1 „r o- 



z l dh*(z) 



D(x;S) 



When we perform a measurement we get an outcome 2° from 
the set D(x;5). The new state of the system z is given by 
formula (fTlj) : 



2 = IE * 2 = 



Tr(x o z°) 



The probability density for the system to be in state z G 
D*(x;5) is 



= (^tliVV, x; S)N(z°, x; S)\ip t ) 



(16) 



where \ipt) is given by J5J and D*(i; 5) is the left translation 
of D(:r; 0") by it (sometimes the notation u*D(x; 8) could also 
be used) defined as 

D^(x;S) = {y G A n |y = u*y°, for some y° G D(x;5)}. 

Expanding (I16|l we get 



E 



x k p\ 



x l p'i 



m=1 J- Pm V Z^m = l x f« 



x J2 -^f^/i^y $k\ PiPj \4>i) 

»>J=1 

_ n -A ir'p°(2°)' _ n Tr(a; op°oz°) 
~ A ^ E„=i z m i& ~ ~A Tr(xopO) 
n Tr(x o p° o (a; -1 * 2)) 



A 



Tr(xop<>) 



(17) 



and if 2 g A n \L>J (x; 5) then fj,(z) = 0. 

We want to emphasize that performing a measurement with 
operators (|15[) on the system at time t gives the new state of 
the system at a later time t + e. For the stochastic process to 
be continuous we want that e — > as 5 — *• 0. 

The easiest way to derive the stochastic differential equa- 
tion for the process X requires the concept of stopping time. 
To give a proper definition we need to consider not just a 
probability space but a filtered probability space. More infor- 
mation on probability spaces, stochastic processes, stochas- 
tic differential equations and Ito's calculus can be found in 

El Q El ESIiai- 

Definition 2 Let (fi, 5, P) be a probability space. A family 
{$(t)}t>o of a -algebras $(t) C is a filtration on the proba- 
bility space (fl, P) if$(s) C ^(t) for s < t. We are going to 
call (Q, $(t)}t>o,$, P) a filtered probability space. 



In this paper we will always assume that the filtered probabil- 
ity space satisfies all the needed hypotheses. 5 is P-complete, 
5t contains all the null sets in J, J t is right- continuous. Also, 
all the stochastic processes are to be adapted to the filtration. 



Definition 3 A process X : R + x Q. -> V (with (V, 2J) a 
measurable space) is adapted if Xt ~ X(t) := X(t,-) is 5t- 
measurable for all t G R+. 

Definition 4 ^4 process X : R+ x f2 — » 1/ (wif/i (V, r) a topo- 
logical space) is continuous ifX(t,io) is a continuous function 
of t for all uj G fl. 

Now we can give the definition of a stopping time: 

Definition 5 A random variable r : fi — ► [0, 00] is called a 
stopping time with respect to {5t}t>o provided 

{t < t} = r _1 ([0,t]) G 5t /or aH i > 0. 

Definition 6 Given a stochastic process X : R+ X f2 — * V 
and a stopping time t the stopped process X T is defined as 

X t — ^min{t,T} 

The stopping times that we will deal with are called hitting 
times. 

Lemma 1 Let X : R+ x fl — > V be an adapted stochastic 
process and S <ZV a measurable set. Then 



t := inf{t > \X t G S} 

is a stopping time. 

The hitting time is just the first time when the stochastic 
process touches the set 5*. 

Now we can derive the stochastic differential equation for 
our stochastic process Xt . Let S G A" and B C dS be a part 
of the boundary dS of 5*. Let Pg (Xo) be the probability that, 
when the stochastic process X with time evolution generated 
by the measurement operators (|15[1 and initial condition Xo 
hits the boundary, dS, it hits it in the subset B rather than 
in dS\B. Let's denote by rgs the hitting time of the process 
of the the boundary dS. 

What we want is to find a second-order differential equation 
for Pg(x). Let's assume that at time t the system state is 
localized at the point x\ Xt = x. Performing the measurement 
with operators (|15[) gives us the state of the system Xt+ £ at 
time t + e which is a random variable with probability density 
distribution given by (|17p . Now we need the following simple 
theorem: 

Theorem 1 Let Xt be the value of the stochastic process de- 
fined above at some time t < tqs- Then if X~ ' is the stochas- 
tic process with time evolution generated by the measurement 
operators M5[) but with initial condition X^* = Xt we have 

pi(x )=p^{x t ). 

The statement of the theorem is trivial because Xj~* — Xt+r- 
From the theorem follows that 

P£(x)=P§(X t+e ). (18) 
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P§ (Xo) is a conditional probability — it's the probability for From this it is a trivial consequence that 
hitting B given that the initial condition is Xo, and thus it 
satisfies the following property 



P§(X ) = / P^{z)¥(Xo 1 (dK{z))) 



Pb {z)@x (dh*(z)) 



= / Pg( z )e Xo (z)dh*(z) 



^fflAffl } 6 aS ^°' (23) 

Now we can take the limit 5 — > of (|21l) . dividing it first 
on both sides by d 2 (<5). For this purpose we first expand the 
integrand in a Taylor series. Because of (|23p , all terms of order 
three or more in the series go to zero as 8 — > 0, and so we need 
to keep only terms of order zero, one and two. To simplify 
(19) the notation we denote p(z) — Pg(z), 6(z) = 9x t+c (z) and 



where Qx a {z) = P(X^" 1 (z)) is the probability distribution of 
Xo and 8x (z) the probability density with respect to the 
Haar measure h*. Substituting this in (1 18 p we have 



P£{x) = P£{X t+e ) = J P£(z)0 Xt+e (z)dh*(z) (20) 

D*(x;5) 

with 9x t+e (z) given by (JTZj) 

nTr(iop°o(x _1 *z)) 

After changing the coordinates z — + a;*2, given that the Haar 
measure is invariant under left translations (h(A) — h(x * A) 
for Vz G A" and A C A n ), d2DJ) becomes 



P§(x) 



P§(x*z)8x t+e (z)dh*(z) (21) 



with 



n Tr(x o p° o z) 
A Tr(xopO) 



(22) 



Equations (|21|l and (|22[) are the starting point for deriving 
the differential equation for Pg(x). We first do an integral 
estimate that is needed to take the limit 8 — > in (|21[). We 
want to estimate the integral 



(x;<5) = y ((* - e) m ) iu "' ik dh*{z) 

D(x;S) 

{z-e) ll ...{z~e) lk dh ir {z). 



D(x;6) 



Using the notation in (|14l) for the i^-th component of z — e 
thought of as a vector in R" we get 



y = x * z. We have 
p(as) =p(y)| e=e 



dp(y) 



dz 



+ 2 l J dz 2 



( z -e) + 0((z-e) 3 ) 



= p(x) + 



dp(y) 



dy 

1 . , dy 



, i dp(v) 



dy_ 

dz 

d 2 P (y) 



2 8y 



+ 0((z-ef). 



dy_ 

dz 



[z-e) 



(z-e) 



Above ■ denotes the usual dot product in R". Substituting 
this in (|21[) while recalling that 



c t+e (z)dh*(z) = 1, 



D(m;6) 



we get 



9y 9 2 p(y) <9y 



+ 2 rr i v az 9y 



2 cty 
+ 0(d 3 (S)A(8)), 



where 



% k {x;S) = J (z- e)® k 6{z)dK{z). 

D(x;S) 



(24) 



(25) 



(z - e) ik < d(x; S) for z G £>(a;; 5) 
and it follows immediately that 



(z - e) ll ...(z - e) 1 " S C(d (5)) as <5 0. 



In the end 



(x;S) < d k (x;S) J dh+iz) 

D(x;S) 

d k (x;S)A(x;S) G G>(d fc (<5)yl(<5)) as <5 -> 0. 



Expanding in a Taylor series. 



6{z) = 6(e) + 



dz 



(z-e) 



= ^ + ^(a!*P )-(»-e), 
and substituting the series into (|25|l . we get 

X fe (x;<5) = -rZ k {x;5) + j(x*p°) ■ Z k+1 (x;S). 



(26) 
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Now we can simplify (|24|l , and after dividing both sides of the 
equation by d 2 (S) we get 



: 



dp(y) 



+ 



dy 
dp{y) 



-Tr I 

2 \ dz 



, 1 dp(y) 



dy 

dz 

dy 

dz 



Zi (x; 8) 
d 2 (8)A(S) 

nZ2(x; 8) 



„ =e d 2 (8)A(S) 



•(x*p°) 



2 dy 
+ <D(d(6)A{6)). 



d 2 P (y) 


dy 


Z 2 (x; 8) 


dy 2 


' dz 

y — x 


z=e <P(S)A(S) 


■Tr(^ 


Z 2 (x;5) \ 


\dz 2 


z=e d 2 (8)A(8)J 



(27) 



The term in the first line of (|27|l is zero because (|13|l can be 
rewritten as 

J (z - e)dK(z) = = Zi(x;8). 

D(x;S) 

Now we can take the limit 8^0. As 

Z 2 (x; 8) 



d 2 (8)A(S) 



€ 0(1) as 8 -> 



it has a finite limit which we denote by rj(x). From its defi- 
nition it follows that this matrix is positive with at least one 
eigenvalue equal to zero, because 



n ~ 

J (z-eY(z~e) k dK(z) 



„ i=l 



(28) 



We will assume that all other eigenvalues of r/(x) are different 
from zero. Because 0(d(S)A(S)) — > as 8 — > 0, we finally 
derive the second-order differential equation for p(x): 



dx dz 



1 f dy 



2„ 



d 2 p dy 
dx 2 dz 



ri(x) ■ (x*p ) 

■V{x] 
,,!,'! ) =0. 

=e 

Differentiating y = x ★ z we get 



2 da; I <9z 2 



<9^ 



and 



»2 i 



dzWz k 



■ ni' (<5] — x 3 ) 



n x (2x J x — x J 8 k — x dj) 



(29) 



(30) 



(31) 



The matrix (|30[) has a pseudo-inverse given by 
and so it satisfies 

£§S i(«f-i) = »-« , )f*f-i 



n 



Substituting these in (|29f) we get the final simple form of our 
equation 



1 a, s d 2 p L . , W/ , <9p „ 



where 



and 



n 2 <9z° 



n a ?(x) ^ 



(x* P y _ P ° 



xi Tr(x o p°) ' 



(32) 



(33) 



(34) 



It's simple to derive the stochastic equations for the process 
Xt having equation (|33p for the probability Pg . We use the 
following theorems to that purpose 

Theorem 2 (Ito's formula for stopping times) Given 
a stochastic process satisfying dX\ = rf(X,t)dt + 
52™ =1 °~] (X, t)dW J , a stopping time r ana! a C 2 ' 1 -function 



ds+ / flu • o- • dW (35) 



u{x,t), it follows that 

T 

u(X t ,t)\l = J (^+Ci 

o u 

where 

^^E^^7 + E4^-E^. 06) 

1,^ — 1 i— 1 fc — 1 

Taking expectations of (|35[) we get 



E(«(jr T ,T))-E(u(jr 0j o)) = E If 



ds 



(37) 



Theorem 3 Given the conditions in the theorem above, the 
C 2 -function u(x) defined on some smooth, bounded domain U 
given by 

u(x) = E(h(X T )) with X = 

satisfies the boundary value problem 

Cu = in U, 
u — h on dll. 

We can now apply this to our case: 

CPS = in S, 



Pg = 1 on dB, 
Pb = on dS\B 



(38) 



with 



Cu = ¥^Sh +h ^ {x) ^- (39) 



dx* 
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From these theorems we see that if we require that the 
stochastic process satisfies the following stochastic differen- 
tial equation it will satisfy equations (|38[) and (|39[) , too 



(40) 



dx' t = y^g iJ '(x t )bj(x t )dt + y^ y a k (x t )dW k 

j = l fc=l 

where o(a;) is the unique square root of g(x): 

n 

ii \ i j 

9 =2^ a h a i- 



(41) 



(|40|) is the equation we consider for governing the evolution of 
the continuous measurement process. We want to point out 
the form of the equation is not accidental. Such equation arise 
when one considers stochastic processes in local coordinates 
on Riemannian manifolds equipped with an adapted (but not 
torsion-free) connection. When one considers Brownian mo- 
tion on manifolds the same equation appears but the term 
in front of dt involves the Levy-Civita connection associated 
with the metric g. For further reading on stochastic processes 
on manifolds one can turn to references [H, QJ, HI EL E2] . 

Now we prove that if the metric g(x) is chosen properly, 
the stochastic process reproduces the strong measurement at 
long times. 

Theorem 4 Let the metric g be invariant under the action 
of A™ on itself and of the symmetric group S n . Then the 
stochastic process satisfying (|40[) starting from Xo = e ends 
in one of the vertices vr^) °f A™ with probability p\ in the 
limit t — > oo. 

The action of S n on A™ is defined as follows — for a G S„ 
and x G A n , 

/ \i cr(i) 

[o~x) = x . 

The components of the vertices vfy are equal to <S|. 
Proof. The conditions for invariance of the metric g together 
with condition (I28p fix the form of the metric very stringently 
— the metric is unique up to a constant conformal factor and 
has the following form: 

n 

g ij (x) = J2 ~ ^W^ty - x 13 ), (42) 

a,0=l 

where the matrix n al3 is constant and equal to the projector 
from R n to A n : 



rf p = (43) 



The stochastic equation (140 f) will have the form 

n 

dx 1 = J2 ^0% - x a )r) al3 x 3 {5 3 p ~ x l3 )b :i (x)dt 

a, /3 = 1 

n 

+ x'(Sa — x a )<j^dW k , 

a = l 

with o is the unique square root of n: 



°~k = Ofc . 

n 



(44) 



(45) 



If we change the coordinates from x to x, defined by 

i = x*p°, (46) 
the equation in these new coordinates looks simple: 

n rt 

dx' = ^i 5 * - x a )o~kdW k = x J2(5'k ~ x k )dW k . (47) 

a=l fc=l 

From this it immediately follows that the stochastic process 
X t in these coordinates is a local martingale; and as it takes 
values in the bounded set A™, it is a martingale. As the 
process is a L p -integrable martingale for any p > 1, it follows 
by the martingale convergence theorem that there exists a 
random variable x^ which is the limit of it as t — > oo. Now 
we will prove that the limit Xoo is localized on the vertices of 
A". Let's denote the m-th moment of x m m (x) — Y2™ =1 (x*) m ■ 
Then using Ito's calculus we get 

n 

dm2(x) = x'(S k — x k )x 1 {8 k — x k )dt 

i,k=l 

n 

+ 2 J2 (x) 2 (Sl~x k )dW k 

i,k = l 

= (tTi2 (i ) — 2m 3 (x) + m 2 (x ) ) dt 

n 

+ (x) 2 (Sl~x k )dW k . 

i,k = l 

Taking the expectation value, we arrive at an ordinary differ- 
ential equation for E(rri2): 



rfE(m 2 ) 
dt 



E(m 2 ) - 2E(m 3 ) +E (m^' 



(48) 



As the limit of it as t — > oo exists, so does the limit of 
E(rri2(ii)). This means that the time derivative of E(rri2(it)) 
goes to zero as time goes to infinity: 



E(m 2 (x oo)) - 2E(m 3 (ioo)) + E (m^ioo)) = 0. 



(49) 



This last equation implies that the range of the random vari- 
able Xoc is a subset of the set of all zeros of the function 
(rri2 — 2m 3 + ml) (x). It is easy to show that the roots of this 
function are exactly the vertices of A" by considering the 
inequality 



(m 2 - 2m 3 + m|) (x) = ^ 



(x*) 2 1-25* + 



> ^(£ l ) 2 (l - 2x l + (x 1 ) 2 ) = ^(£ l ) 2 (l - x) 2 > (50) 

i = l i=l 

with equality when x G {^(i), .. .,«(„)}. This proves that the 
stochastic process ends at one of the vertices of A" when time 
goes to infinity. 

Let's denote the probability distribution function of Xoo by 
m(x). Then 



m(x) 



q k , for x = v (k) 

0, for x ^ H( fe ) for Vfc : 



(51) 



for some q = (qi, ... , q n ) G A™. As x t is a martingale it follows 
that 



E(i ) =E(xt) =E(5 00 ). 



(52) 
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In the ^-coordinates the process starts from xo = e. By (|46 
the initial condition in the ^-coordinates is xo = p°. Then 



p° = E(x ) = E(x^) = 1kv\ k) = ^2 q k Sl = q t . (53) 

fc = l k = l 

This proves the theorem. 

Now we can describe the evolution of the state of the system 
when it is subjected to our feedback control scheme. The 
equations are a special case of equations (|61[1 with Mj = Pj\ 

3=1 

1 n 

+ ^J2(P-{h))\4>t}dW\ (54) 



dx l = g 13 ( ) dt + a a dW° 



(55) 



with ( Pj ■) — {ipt \ Pj \ ipt) ■ The notable thing in these equa- 



tions, that distinguishes them from the equations describing 
the the decomposition of a generalized measurement, is that 
feedback is not needed. The state of the system \ipt) evolves 
in accordance with equation (|54l) which has no dependance on 
the vector x. As we will see below, that is not true for gener- 
alized measurement. In that case the equations for \ip t ) and 
x are coupled in a nontrivial way, and feedback is necessary. 



IV. CONTINUOUS PROCESSES FOR 
GENERALIZED MEASUREMENTS 

Building from the above process on the classical state space 
A", it is now easy to prove the existence of a stochastic pro- 
cess on the quantum state space that decomposes any kind 
of generalized measurement, given by measurement operators 
Mj (j — l,...,n) satisfying the usual completeness relation 
S?=i MjMj = /, into continuous measurements. (It is not 
necessary that these measurement operators commute.) To 
that purpose we will use the fact that the homomorphisms of 
stochastic differential equations(the transformations that pre- 
serve the structure of the equations) must be at least twice 
differentiable maps. We will give an example of a map M 
taking points in A n and mapping them to operators acting 
on the state space of our quantum system. This map should 
satisfy the following three requirements: it should be twice 
differentiable in A n , equal to the identity I at the identity e 
of A n , and equal to the measurement operators at the vertices 
of A n . 

Let Mj = UjLj be the left polar decomposition of Mj, so 
Lj are positive and Uj are unitary. The map M will be the 
product of two maps T(x) and A(a;): the first involving the 
unitaries Ui, and the second the positive operators Lj. If the 
hamiltonians corresponding to the unitaries Uj are Hj, then 
we can choose the map T(x) to be 



T(a;) = exp 



2=1 



(56) 



We can readily see that, as constructed, the map has all the 
required properties, and also has the property that T(x) is 



unitary for all x G A™. A(x) is given by 
k{x) = /*(*) (fVLjL. 



3=1 



(57) 



with f(x) = l + nEL^fl 



The square root above 



is well-defined because the the expression in parentheses is 
a positive operator for every x £ A". We also note that 
A(i) is invertible for every x G A™. This follows from the 
completeness property satisfied by the {Lj}: X^=i A7 = I- 
The map Mix) is given by T(a;)A(a;). 

The quantum state p will respectively evolve according to 
the following formula 



Pt = 



M{xt)poM\x t ) 
Tr(M1{x t )M(x t )po) 



(58) 



It is obvious that the constructed map is far from unique — 
there are infinitely many maps that satisfy the three required 
conditions. This shows that there are numerous ways to de- 
compose a strong generalized measurement into weak mea- 
surements of the type that we are considering. 

In the very special case where the measurement operators 
are positive and commuting, we will now show that the state 
evolves according to a generalized quantum state diffusion 
equation. As the measurement operators are positive, the 
unitaries in the their polar decomposition are all equal to the 
identity operator. Thus the map M(x) is just A(a;). A pure 
state will evolve by the analog of (f58|) : 



A(x t ) |Vo) 



{ipo\ A 2 (xt)\ipo) 2 



(59) 



As the measurement operators commute, it is easy to take 
derivatives, and then use the Ito rule to get the quantum dif- 
fusion equation for \ijjt) ■ In the following we use the Einstein 
summation rule over repeated upper and lower indices: 



d\i>t) 

d\A) 
dx l 

dxWx k 



dx i 

J 

2 



+ 2dxWx<< a ' 



M 2 



M? 



x m M%, 
( 



x m Ml 



M£Mj 

(x m Miy 



+ 



Mj{M£ 



x m M2 ( x l M? 



+ 



Mf) Mi 



Mj 



x m Ml.) x l Mf 







x m M%, 



Ml 



dx' — g XJ bjdt + a a dW a 



x m Ml 



-dt + a a dW° 



dx 3 dx k = g 3k dt, 
with (Mf) = (ip \ M] |Vo) = (Vol M\Mj |Vo) • It's easy to 





see that , ^ Mj .\\ 
Denote 



A 1 



where (Mf) = (Vt| Mf \ip t ) . 



(60) 



/ \ 3 / ' 

Mf 



x m Ml 
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Putting this all together, we arrive at the following coupled 
stochastic differential equations: 

d\ip t ) = -lg jk (Aj - {A 3 ))(A k - (i fc )) dt 

dx l = g l] (Aj\ dt + a\dW a . (61) 
From here we can easily derive an equation for pt = \ipt) {ipt\'- 

dpt = g 3k [QjPtQk - \{QkQ.uPt}\ dt + {Q 3 , Pt }aUW a , 

dx l = g lJ Tr (Ajpt) dt + a l a dW a , (62) 
where 

Qi(pt,x) = \ (Aj{x) - Tr (ij(x)pt)) . (63) 

V. CONCLUSIONS 

In this paper we have constructed stochastic processes that 
are continuous decompositions of a discrete, instantaneous 
measurement. The state of the system evolves in accordance 
with stochastic differential equations, which in the case of 
generalized measurement involve feedback based on the mea- 
surement history. We have proven that in the long-time limit 
this process arrives at the same final states with the same 



probabilities as the strong measurement. This gives us the 
ability to think in terms of this continuous process when we 
consider measurements, and a natural way to express state- 
ments involving measurements in the language of stochastic 
calculus. For example, we can easily differentiate with re- 
spect to the process using Ito's calculus. One application of 
this idea has already been made to entanglement monotones, 
giving new differential conditions for them, as shown in [3J. 

A related and still unsolved question concerns the converse 
problem to the one considered in this paper: namely, if we 
are able to perform a certain restricted class of operations 
or if we have control over certain parameters of a quantum 
system, to what extent we can utilize this freedom to set up 
an experiment in which the system's evolution is governed by 
the model considered above? Given some class of weak mea- 
surements that can be experimentally performed, what class 
of generalized measurements can we generate? This question 
has relevance when one starts exploring an experimental real- 
ization of the proposed type of continuous measurement, and 
is the subject of ongoing research. 
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